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The question about the existence of a structural glass transition in two dimensions is studied 
using mode coupling theory (MCT). We determine the explicit d-dependence of the memory func- 
tional of mode coupling for one-component systems. Applied to two dimensions we solve the MCT 
equations numerically for monodisperse hard discs. A dynamic glass transition is found at a critical 
packing fraction v^f^^ = 0.697 which is above <Pc^^ — 0.516 by about 35%. (/j^ scales approxi- 
mately with ipr(.p the value for random close packing, at least for d—2, 3. Quantities characterizing 
the local, cooperative 'cage motion' do not differ much for d = 2 and d = 3, and we e. g. find 
the Lindemann criterion for the localization length at the glass transition. The final relaxation 
obeys the superposition principle, collapsing remarkably well onto a Kohlrausch law. The d = 2 
MCT results are in qualitative agreement with existing results from MC and MD simulations. The 
mean squared displacements measured experimentally for a quasi-two-dimensional binary system of 
dipolar hard spheres can be described satisfactorily by MCT for monodisperse hard discs over four 
decades in time provided the experimental control parameter F (which measures the strength of 
dipolar interactions) and the packing fraction ip are properly related to each other. 

PACS numbers: 64.70.Pf, 61.20.Lc, 61.43.Fs 



I. INTRODUCTION 

Static and dynamic behavior of macroscopic systems 
depends sensitively on the spatial dimension d. For ex- 
ample, one-dimensional systems with short-range inter- 
actions do not exhibit an equilibrium phase transition. 
In two dimensions there is no long range order if the 
ground state exhibits a spontaneously broken continu- 
ous symmetry, and Anderson localization occurs for al- 
most all eigenstates of a disordered system for d—\ and 
d=2, but not in (i=3, if the disorder is small. Critical 
exponents at continuous phase transitions depend on di- 
mensionality. Concerning dynamical features it is known, 
for instance, that the velocity autocorrelation function of 
a liquid exhibits a long-time tail proportional to t^'^/'^. 
Consequently, the diffusion constant is infinite for d <2. 
These few examples demonstrate the high sensitivity of 
various physical properties on the dimension d. 

Let us consider a liquid in d—3. If crystallization 
can be by-passed a liquid undergoes a structural glass 
transition. Although not all features of this transition 
are completely understood, recently significant progress 
has been made concerning its microscopic understanding. 
Following many decades of several phenomenological de- 
scriptions with less predictive power, the mode coupling 
approach introduced in 1984 by Bengtzelius, Gotze and 
Sjolander , has led to a microscopic theory of the struc- 
tural glass transition. 
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This theory, called mode coupling theory (MCT), has 
been discussed theoretically in great detail by Gotze and 
his coworkers (see Ref. [2] for a review). Its numer- 
ous predictions were largely successfully checked by ex- 
periments and simulations [3j |4j. The main prediction 
of MCT is the existence of a dynamical glass transition 
at which the dynamics changes from ergodic to noner- 
godic behavior. Thermodynamic (equilibrium) quanti- 
ties, e.g. the isothermal compressibility and structural 
ones like the static structure factor S{q), do not be- 
come singular at the glass transition singularity of MCT. 
Hence, the MCT glass transition is of pure dynamical 
nature. It can be smeared-out by additional relaxation 
channels, and then marks a crossover [21 [3] ■ 

A microscopic theory predicting a structural glass 
transition with pure thermodynamic origin was derived 
by Mezard and Parisi f5|. Their replica theory is a first 
principles approach which yields a so-called Kauzmann 
temperature Tr- at which the configurational entropy 
vanishes. Tk is below Tc, the MCT glass transition tem- 
perature. In two dimensions, Tk may mark a crossover 
[Bj. For a review of both microscopic theories as well 
as phenomenological approaches to the structural glass 
transition the reader may consult Ref. [7|. 

An important question is now: "What is the depen- 
dence of the structural glass transition on the spatial di- 
mensionality" ? This question has already been asked by 
several people some time ago. Before we come to a short 
review of this work, let us consider monodisperse hard 
spheres and hard discs in d=3 and d=2, respectively. 

The most dense packing of four hard spheres cor- 
responds to a regular tetrahedron. However, three- 
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dimensional space cannot be covered completely by regu- 
lar tetrahedra, without overlapping. This kind of geomet- 
rical frustration is absent in two dimensions. The densely 
packed configuration of three hard discs corresponds to 
an equilateral triangle. Since the two-dimensional plane 
can be tiled completely without overlap by equilateral 
triangles, there is no frustration. Therefore one may be 
tempted to conclude that there is no structural glass tran- 
sition in two dimensions. However, the link between frus- 
tration and glass transition has proven subtle. Experi- 
ments [8J and simulations ^ [TD] of monodisperse hard 
spheres in three dimensions have shown crystallization. 
Hence, the existence of geometric frustration is not suffi- 
cient for glass formation. What one needs is bi-dispersity 
or polydispersity. 

Santen and Krauth have performed a MC simulation 
of a two-dimensional system of poly disperse hard discs. 
The polydispersity has been quantified by a parameter 
gPoi_ 'pj-^gjj. results clearly demonstrate (i) the absence of 
a thermodynamic glass transition and (ii) the existence 
of a dynamic glass transition at a critical packing fraction 
<Pc(e''°')- The kinetic glass transition shifted outside the 
region of crystallization for > eP°|^ « 10% [TTKH]. 
Furthermore, the diffusivity was found to be consistent 
with the MCT result [21 [3]: 

where 7 w 2.4 and ip™^ w0.80. This critical value agrees 
with what has been found for a related system by Doliwa 
and Heuer (see Figure 2 in Ref. fW). The absence of 
a thermodynamic glass transition has been strengthened 
recently by Donev et al. |14j for a binary hard-disc mix- 
ture. 

There are a few investigations of glass formation in 
two-dimensional systems with soft potentials. Langon 
and Chaudhari [15J studied a binary system with mod- 
ified Johnson potential. They found that the structural 
relaxation time seems to diverge when approaching a 
critical temperature. Similar behavior was observed by 
Ranganathan |16j for a monodisperse Lennard-Jones sys- 
tems and by Perera and Harrowell '17] for a binary mix- 
ture of soft discs with a 1/r^^ potential. It is surprising 
that the intermediate self-scattering function S'-^''{q,t) of 
the monodisperse system |16j . exhibits strong stretching, 
one of the characteristics of glassy dynamics. However, 
5^**^ (q, t) does not produce a well pronounced plateau [TB] 
under an increase of the density, i.e. the cage effect does 
not become strong enough. This is quite different to the 
binary system |17j . S^'^\q,t) develops a two-step relax- 
ation process upon supercooling with a well pronounced 
plateau over 4-5 decades in time, at lower temperatures. 
This behavior is qualitatively identical to that found for, 
e.g. the Lennard-Jones mixture investigated and ana- 
lyzed in the framework of MCT by Kob and Anderson 
[H]. The authors of Refs. [HI [HI [HI El HSl HZ] 
conclude that there is a structural glass transition in two 
dimensions. Their conclusion is supported by recent ex- 
periments on colloidal particles with repulsive dipolar 



interactions in two dimensions [TO]. Since these simu- 
lational and experimental findings strongly resemble the 
MCT predictions obtained for d=?>, it is important to 
apply MCT to two-dimensional liquids. MCT has been 
applied to the two-dimensional Lorentz model of over- 
lapping hard discs [20 and a charged Bose gas with 
quenched disorder and logarithmic interactions at zero 
temperature [21] , but to our best knowledge not to a 
two-dimensional liquid-glass problem with self-generated 
disorder. To accomplish this is the main motivation of 
the present contribution. 

The outline of our paper is as follows. The MCT equa- 
tions for arbitrary dimensions and the major predictions 
of MCT will be presented in Sec. II. The theory re- 
quires as only input the static structure factor, which is 
computed in Sec. III. In Sec. IV we apply MCT to a 
two-dimensional system of monodisperse hard discs and 
will demonstrate that there is a dynamic glass transition. 
The dynamic behavior close to that transition is qualita- 
tively identical to that of monodisperse hard spheres in 
three dimensions. It will also be shown that the MCT- 
result for the time dependent mean squared displacement 
describes the experimental result of Ref. [H] for both 
species rather satisfactorily over 4 decades in time. The 
final section V contains a short summary and some con- 
clusions. 



II. MODE COUPLING EQUATIONS 

In this section we will shortly review the MCT equa- 
tions for the collective and tagged particle correlator 
of density fluctuations of a one-component liquid and 
will present the properties of their solution close to 
the glass transition singularity. The only dependence 
on dimension d comes through the integrations element 
(27r)~'^(i'^fc ~ {2n)~'^¥^~'^dk which appears in the mem- 
ory kernels. 

MCT provides equations of motion for the normalized 
intermediate scattering function 4'q{t) and the tagged 

particle correlator (j3^q\t). The mathematical structure 
of these equations does not depend on d. For Brownian 
dynamics which is appropriate for colloidal systems they 
read for d arbitrary 

t 

7,</),(t) + (t)g{t) + J dt'niqit - t')^q{t') = (1) 


with the memory kernel mq{t) containing fluctuating 
stresses and playing the role of a generalized friction co- 
efficient. It arises because the density fluctuations cap- 
tured in 4>q{t) are affected by all other modes in the sys- 
tem. In MCT, one assumes that the dominating contri- 
butions at long times are given by density pair fluctua- 
tions and approximates (in the thermodynamic limit) 

mqit) = Tq[Mt)] = F(q,k,p)0fe(t)^p(t) (2) 
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The vertices express the overlap of fluctuating stresses 
with the pair density modes, and are uniquely determined 
by the equilibrium structure 

V^(q,k,p)= ^;?i^[q.kcfc + q.pcp]2<5(q-k-p) (3) 
z q 

where n is the number of particles per d-dimensional vol- 
ume, Sq the static structure factor and Cq the direct cor- 
relation function related to Sq by the Ornstein-Zernike 
equation. 7^ is a characteristic microscopic time scale. 
The reader should note that the vertices, Eq. (|3|, have 
been approximated by neglecting static three-point cor- 
relations. 

The corresponding equations for the tagged particle 
correlator and d arbitrary are of the same form 

t 

-fi^^^^Ht) + 4'\t) + J dt'm^^\t - t'),p<q'\t') = (4) 


with 



and 



W^)(q,k,p) = n5fcg-4(q.k)2(4^Vj(q_k-p) . (6) 



c\^' = {pq * Pq) /{nSq) IS the tagged particle direct corre- 
lation function. If the tagged particle is one of the liquid's 

particles it is c^^"^ = Cfe. 

The static correlation functions S^, Sp, Sq,Ck,Cp and 

the correlators 4>k{t), 4>k\'^)^ 4>p{t), 4>p\t) depend on 
q = |q|, k = |k| and p = |p|, only, due to the isotropy of 
the liquid and glass phase. Therefore the d-dimensional 
integrals in Eq. Q and Eq. ([s]) can be reduced to a two- 
fold integral over k and p. This will make explicit the 
d-dependence of the vertices. As a result one obtains, 
similarly to the case in d = 3 [I], 
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iq' + e-p')ci'^] Mt)4'\t) (8) 
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with which is consistent that the "one-dimensional unit- 

^ sphere" is an interval of length two with a "surface" con- 

Q _ 2'""^ /Q\ sisting of two points. 

^2' The behavior for q ^ of both functionals can be 

k+q 

the well-known result for the surface of a d-dimensional obtained by a Taylor expansion of / dp{---). Although 
unit-sphere. T{x) is the Gamma function. Note that k-q 
Eq. (Igj) yields (with r(i) = tt^/^) straightforward it is rather tedious. Alternatively, one 

can start directly from Eqs. ([2]) and ([s]), in order to 
f^i = 2 (10) obtain for q^Q 
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(11) 
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and 

oo 




Note that ^o[0fe(^)] exists, whereas Tq'^\4'k{t), (j)''^\t)] di- 
verges like q~^. This divergence is related to the absence 
of momentum conservation for the tagged particle. 

Taking d=3 in Eqs. (0, ([8]), ([TT]) and ^) one arrives 
at the well-known representations of the memory kernel 
for finite g and g ^ [TJ [221 [23]. Note that Refs. 
and US] already present the integrals in Eqs. ([7| and ^ 
in discretized form. 

With the knowledge of the number density, of the 

(s) 

static correlators 5*^, Cq and Cq as function of the ther- 

(s) 

modynamic variables and of 7g and jq one can solve 
Eq. ([T]) and Eq. Q for initial conditions 4>q{0) = 1 and 

(j)[^\o) — 1. There exist several quantities which char- 
acterize the solutions. These quantities can be found in 
Refs. [2], [22] and [23] • In order to keep our presentation 
self-contained as much as reasonable we discuss those for 
which results will be reported in the next section. We 
start with the glass transition singularity. At the glass 
transition the nonergodicity parameters (NEP) 



fq = lim 



(13a) 



change discontinuously from zero to a positive nonzero 
value, smaller or equal to one. The corresponding quan- 
tity 



/(^)- lim 4^Ht) 



and the so-called exponent parameter, A, which obeys 
< A < 1. The second term on the r.h.s. of Eq. (15) is 
the leading asymptotic result for fq, and Kq + n yields the 
next-to-leading order correction. C, A, and Kq + k fol- 
low from J^q[fk] and its derivatives with respect to fk at 
tpc- Particularly, is the right eigenvector Cq belonging 
to the largest eigenvalue -E'inax(<p) of the stability matrix 
{dTq[fk\/dfk) at the critical point. i?max(</') is not de- 
generate since the stability matrix is non-negative and 
irreducible. At the critical point the maximum eigen- 
value becomes one, i.e. (fc can be determined from the 
condition -B,nax('rf3c) = 1. 

At the critical point (t>q{t) decays to the plateau value 

fc 
J q 1 



< /q < 1. Its time dependence is given by 



4'qit) = f^ + hqit/to)-''{l + [Kq + Kia)]it/tor''} , t » io ■ 

(18) 

Kq (not to be confused with Kq) and k(x) are again 
determined by J^q[fk] and its derivatives at (p — (pc- They 
are a measure of the next-to-leading order contribution 
with respect to the leading asymptotic result 



<l^qit) = fq 



hqit/toY 



(19) 



(13b) 



the critical law. This critical decay occurs on a time 
scale t much larger than a typical microscopic time to- 
The exponent a is determined by A, only. 

For e < two cr-dependent, divergent time scales exist 



can change discontinuously at the same point or in a con- 
tinuous fashion at higher densities or lower temperatures. 
Both NEP fulfill nonlinear algebraic equations ^ 

Y^=^q[fk] , = ^!,'Hhjt^] . (14) 

^ Jq 1 - /<; 

and /q*'''^ are the NEP at the critical point, e.g. at 
n = rii-. Since we will apply MCT to d-dimensional hard 
spheres with diameter D we use in the following the pack- 
ing fraction ip = n^ld-i{D/2)'^ /d. Above, but close to 
ip = (pc i.e. for < e = {tp — ipc)/^c <C 1 it is: 



fq = fq + hqW<j/{l - A) + <j{Kq + k) / Vl^] (15) 

with the critical amplitude 

K = i^- fqf^q , (16) 

the separation parameter 

a{e) = Ce + ©(e^) , (17) 



h \e\- 



^0 



and 



f' - f' 
— t'fi 



e < 



(20a) 



(20b) 



with 7 = ^ + 55, and = to/B^^'', where i? is a con- 
stant. The so-called von Schweidler exponent, b, follows 
from A, only. (t>q{t) exhibits a two-step relaxation. The 
relaxation for t/tcr <C 1 to the critical plateau value 
follows from Eq. ( [T8| by replacing (t/to) through {t/t„) 
and the decay from that plateau to zero is initiated by 
the von Schweidler law for t„ <^t <^t'^ 

m - fq-hq(t/0\^-[Kq + ^^{-b)]{t/t',f} . (21) 

Kq + n{—b) determine again the next-to-leading order 
contribution. 

For e > there is a single relaxation process, only. 
4>q{t) relaxes for t/tc <C 1 like for e < 0, and finally the 
plateau value fq is reached by an exponential long time 
decay. 
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For e < 0, the final or a-relaxation process describes 
the decay of the correlators from the plateau down 
to zero. Asymptotically close to the transition, the func- 
tional form of the a-process is given by a master function 
(j>q{t) of the rescaled time t = t/t'^ via 

= + e ^[^Hi) + 0{e^) for e -> 0- (22) 

The master function (j)q obeys an equation similar to 
Eq. ([ij with vertices evaluated right at the critical point, 
e — 0. Thus it does not depend on separation e and 
control parameters, and Eq. (22 1 expresses the often ob- 



served '(time-temperature) superposition principle' [2J[3]. 
The von Schweidler series Eq. (211 gives the short time 
behavior of (f>q for t — > 0, and the corresponding result 

0) ^ hqBiCi-\ with C 



for the correction is: (t>'^\i 



from Eq. (17 1 and Bi a known constant. 



Similar leading order and next-to-leading order con- 
tributions can be derived for the tagged particle corre- 
lator (j)^q\t), and e.g. the mean squared displacement 
5r'^{t) = ((r(<) - r(0))2) Since 



2d, 



(23) 



the long- wave limit of Eq. ^ yields after integration with 
respect to t: 

5r^{t) + D['n dt'm'^o'\t-t')Sr'^{t') = 2dDi'h (24) 

where the memory kernel mQ^\t) follows from Eqs. (|5 
and Eq. 



mi'\t) ^ limg2^(^)[0,(t),<^i^)(i)] 

oo 



d{2TxY 



dkk''+'Sk{ci''>fMt)<l>'k'{t){25) 



is) I 



Furthermore we have used: 



(26) 



(s) 

with Dq the short-time diffusion constant of the tagged 



particle. Eqs. (13b I and Eq. (23 1 imply that the long- 
time limit of 5t {t) 



lim 6r^{t) 

t — *-CXD 



2dri 



(27) 



is related to the tagged particle's localization length 
given by: 



lim 

q^Q 



is) A I 

In the liquid phase where fq — 0, Eq. (27 1 gives — 



(28) 



oo, i.e. the particle is delocalized. While becomes finite 
at the glass transition. 



Besides the correlators (f>q{t) and (f>q'\t) one can 
also study the corresponding susceptibilities Xqi'^) ^tnd 
Xq'^\i^), respectively. Similar asymptotic laws and next- 
to-leading order correction exist for them [22l [23] . Inde- 
pendent on whether the correlators or their susceptibil- 
ities are considered the dependence on d of the leading 
and next-to-leading order terms enters only via the d- 
dependence of J^, (Eq. ^) and J^^'^ (Eq. 



III. STATIC STRUCTURE 

In this section we consider the calculation of accurate 
equilibrium structural correlation functions for the hard 
disc system. Within MCT, all information regarding the 
interparticle interactions is contained in the static struc- 
ture factor which enters the memory function vertices 
Eqs. (|3| and the interaction potential does not enter 
explicitly in the MCT equations. Experience with MCT 
calculations in three dimensional systems has shown that 
the location of the glass transition somewhat depends on 
the details of the input structure factor, particularly the 
height of the main peak. Different approximate theories 
for the static structure lead to varying values for e. g. the 
critical packing fraction 4>i^~^'' [H]. We have therefore 
considered a number of approximation schemes for the 
two dimensional static structure factor in order to obtain 
the best possible values for the description of the ideal 
glass transition. The quality of the various approxima- 
tion schemes is assessed by comparison with computer 
simulation data. (Monte Carlo simulations of 6 10"* par- 
ticles, as well as event driven molecular dynamics simula- 
tions of 1089 particles were performed, both in the NVT 
ensemble [26].) 

Integral equation theories based on the Ornstein- 
Zernike (OZ) equation provide a powerful method to cal- 
culate the pair correlation functions for a given interac- 
tion potential [25] . The OZ equation is given by 



h{r) = c(r) + n d'^r' c(r') h{r - r'), 



(29) 



where h{r) = g{r) — 1. This expression must be supple- 
mented by an additional (generally approximate) closure 
relation between c(r) and h(r). For hard spheres in d- 
dimensions the most widely used closure is the Percus- 
Yevick relation g{r < 1) = 0, c{r > 1) = 0. In odd 
dimensions the resulting integral equation can be solved 
analytically for the direct correlation function c(r). In 
even dimensions there exists no analytic solution and full 
numerical solution is required [27j . Efforts have been 
made to approximate the numerical PY data by analytic 
forms [2H1 [in] but in all cases these fail to reproduce ac- 
curately the detailed structure of the numerical solution 
at high densities. The formally exact closure to the OZ 
equation for systems with pairwise interactions is given 

by 

h{r) = -1 + exp(-/3M(r) -I- h{r) - c(r) -I- B(r)), (30) 
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where u{r) is the pair potential and B{r) is the bridge 
function, an intractable function representing the sum of 
the most highly connected diagrams in the virial expan- 
sion. Setting B(r) — recovers the familiar hypernetted- 
chain approximation. 




FIG. 1: Comparison between theoretical MHNC structure 
factors and simulation for hard discs at packing fractions (p = 
0.5 (+), 0.628 (o), and 0.68 (A); q is given in units of 1/L>, the 
inverse diameter. Inset (a) concentrates on the vicinity of the 
main peak and gives additional comparison with Baus-Colot 
and PY theories for packing fraction <j) = 0.628. Inset (b) 
demonstrates how the MHNC theory correctly captures the 
asymptotic behaviour for large q values for packing fraction 
6 = 0.628. 



The modified-hypernetted-chain (MHNC) approxima- 
tion is to take B{r) from the PY theory solved at some 
effective density n* , different from the true system den- 
sity n, and to treat this as a variational parameter to en- 
sure thermodynamic consistency between the virial and 
compressibility routes to the pressure. A detailed de- 
scription of the MHNC equation can be found in j30J. 
The steps taken in solving the MHNC equation are, (i) 
numerically solve the PY equation at density n* , (ii) use 
Eq.(|30| to find B{r) = Bpyir; n*), (iii) solve Eqs.( p9p0l ) 



with this bridge function, (iv) calculate the pressure from 
the virial and compressibility equations, (v) adjust n* 
until the two pressures are equal. In three dimensions 
it is generally recognized that the MHNC approximation 
provides a highly accurate description of the pair corre- 
lations for the hard sphere fluid, significantly improving 
upon the PY theory. We find that the same is true in the 
case of two dimensional hard discs. At low densities the 
MHNC Sq lies very close to the PY result. As the density 
increases discrepancies begin to arise, particularly in the 
region of the main peak, with the MHNC in closer agree- 
ment with simulation. Both the MHNC and PY theories 
are significantly more accurate than the analytical Baus- 
Colot expression [28] . Figure [l] shows a comparison be- 
tween the MHNC Sq and the simulation results. The level 
of agreement is very satisfactory and the MHNC shows 



clear improvement over the other theories investigated. 
The only notable deviation from the simulation results 
is the height and width of the second (and third) peak, 
which is overestimated (respectively underestimated) by 
the MHNC theory. It is known that upon approaching 
the crystallization phase boundary (located at (pp = 0.69 
for hard discs) a shoulder develops on the second peak of 
the structure factor; a feature which has been interpreted 
as an indicator of approaching crystallization The 
development of the shoulder also suppresses the height of 
the second peak to some extent and leads to a small shift 
in the location of the third peak. The MHNC theory, like 
the PY and all other standard integral equation theories 
do not contain information about crystallization and thus 
predict fiuid-like structure at, and beyond, the freezing 
transition. While this property leads to some discrepancy 
with simulation results at high density, it makes such the- 
ories ideal for calculating the fluid-like structure factors 
required as input to the MCT, where we assume crystal- 
lization to have been suppressed. We can thus proceed 
with confldence using MHNC structure factors as input 
to the MCT. 



IV. RESULTS 

In this chapter we will apply MCT to a two- 
dimensional system of monodisperse hard discs with di- 
ameter D. We will solve the MCT equations ([T]), Q and 
( 24 1 and will present results for those quantities discussed 
in the last section. As input we will use the static struc- 



|T} and 



we take 



ture factors Sq obtained from the MHNC approach. This 
result is presented in Figure [2] for three different packing 
fractions close to the glass transition and compared with 
the corresponding Percus-Yevick result for hard spheres 
in d=3. For instance, for (p = (p^^'^ and (p — fi^'^, i-e. for 
e = 0, the peaks are more pronounced in d = 2. In par- 
ticular the main peak is more narrow and higher for d — 2 
than for d = 3. The direct correlation function Cq follows 
from the Ornstein-Zernike equation. Because we choose 
the tagged particle as one of the liquid particles we have 

(s) _ 
Cq — Cq. 

For the friction coefficients in Eqs. 
7, = Sq/{D^o'\^) and 7^'^ = l/{Di''q') (cf. Eq. (|26|)) 
and choose as our unit of time t^^^^ = Dq^^ /{lOD^); 
in the following all times will be given as rescaled ones, 
^I^(BD) ^ For the numerical solution of MCT equation 
one has to discretize q. A compromise between a fine grid 
and computation time is required, as the computations 
scale with number of gridpoints . We choose a grid 
with M — 250 gridpoints and a high q cutoff of 50/-D. A 
higher cutoff has only a very small effect on the critcal 
packing fraction (< 10^"*). The effect of the number 
of gridpoints is more sensitive due to the form of the 
Jacobian of the transformation to bipolar coordinates. 
In this paper the integration is substituted by a rule that 
can be called a modified trapezoid rule. The value of 
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the function to be integrated is not taken in the middle 
of the interval [nh, {n + l)/i] but at (n + 0.303)/i. By 
using this rule on gets the best discrete description of 
the Jacobian hence it is used here. The difference in </? 
between M = 500 and M = 250 is then < 10"^ leading to 
a system close enough to continuum. Then the solution 
of Eqs. (jlj, Q and (24 1 will be performed by use of a 
decimation technique [32J . 





FIG. 3: Non ergodicity parameter of coherent (2d diamond, 
3d dashed line) and incoherent (2d circles, 3d solid line) cor- 
relators at critical packing fraction tpl'^ = 0.697 and ip'i'^ = 
0.516. The q-values for q < 1 are not included since they can 
not be determined accurately for numerical reasons. The q=0 



results are from the analytic expansions Eqs. (Ill and (12 1 
and are included as triangles. 



FIG. 2: (Color online) Comparison of 2d MHNC structure 
factors (black/dark) to 3d PY (magenta/light). The pack- 
ing fractions correspond to e = (solid lines), e — —10"*''^ 
(dotted lines) and e — -1-10"*'''^ (dashed lines). The critical 
packing fractions are 0.697 for 2d and 0.516 for 3d. 



The search for the glass transition singularity can be 
done either by an iterative solution of the non-linear 
equations (14) or by calculation of i?niax(v')- In this pa- 



per a simple bracketing algorithm is used starting from 
two points were point A yields a finite NEP for q near the 
peak position and point B has fq = 0. The next point C 
is taken in the middle of the interval. If it yields finite fq 
the next point is taken between A and C, if = the 
other interval [C,B] is taken. This procedure is continued 
until the critical packing fraction tpc is determined to a 
precision of 10^^". As a result we have found 



ipi=^ ^ 0.696810890(317) 



which is above the value for d — S 



=3 



0.51591213(1) 



(The denoted accuracy will be required to reliably com- 
pute e in the following.) Similar to the three-dimensional 
system, the collective and self part of the density fluctu- 
ations become nonergodic at the same critical packing 
fraction lysf^^. The corresponding critical nonergodicity 
parameters are shown in Figure |3] 



While almost no difference between incoherent NEPs 
for d = 2 and d — i can be observed, more pronounced 
maxima appear in the coherent NEP at higher wavevec- 
tors for the lower dimension. Regions of rather abrupt 
q-dependences in fq should be observable experimentally. 
Two length scales appear to be involved in fq. While the 
average particle distance, connected to the main peak in 
Sq, somewhat differs from d = 3 to rf = 2, the localisation 
length, which dominates the incoherent NEP, is insensi- 
tive to dimensionality. The change of fq when stepping 
down in dimension thus can not simply be scaled away. 

An important observation in the numerical solution of 



Eq. (14) concerns the convergence of the required inte- 
grals. We find here, and for all other integrations per- 
formed, that convergence at small and large wavevectors 
holds. We interpret this as indication that the MCT glass 
transition in d = 2 describes a local phenomenon not af- 
fected by long range correlations, which might sensitively 
depend on dimensionality. 

All wavevector dependent structure functions describ- 
ing the glassy structure and its relaxation ('cage effect') 
are summarized in Fig. [4] The critical nonergodicity pa- 
rameter and the critical amplitude hq is depicted in 
Figure 4b and t he n e xt-to -leadi ng order amplitudes Kq 
and Kq (cf. Eqs. ([l5l, ([is]) and (|2l])) in Figure|4[;. The in- 
dependence of the shown quantities follows generally the 
one of Sq . Only hq exhibits the opposite variation around 
the main structure factor peak. Figure |4^ specifies three 
typical wave numbers for which results will be discussed 
below. Comparison with the d — "i result (Figure 2 in 
Ref. [22]) reveals qualitative similar g-dependence of 
hq, Kq and Kq. Like for Sq, the (/-variation of all quan- 
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titles is more pronounced in d = 2 than in d = 3. 
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FIG. 4: (a) Structure factor Sq as function of wave vector q 
for = ~ 0.697 (solid), ip ^ 0.729 (dotted) and (p « 0.664 
(dashed); the latter correspond to e = ±10"'*'''^. The arrows 
mark the wave vectors qi = 6.46, 52 = 10.06 and — 18.26. 
(b) Critical NEP (diamonds) and critical amplitude hq 
(squares), (c) The amplitudes Kg (triangles) and Kq (circles). 
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FIG. 5: Non ergodicity parameter fq for qi = 6.46, q2 = 10.06 
and ga = 18.26 (solid lines). The leading asymptotes (dashed) 
describe fq — fq within 10% up to the values marked by di- 
amonds (0.005, 0.075, 0.0015). The next to leading asymp- 
totes (dotted) are within 10% for e up to the values marked 
by circles (0.03, 0.033, 0.05). The values for the correction 
amplitudes are C = 2.08, a = 1.18, hi = 0.337, /i2 = 0.654, 
hs = 0.508, Ki = -1.86, K2 = -0.456 and K3 = 0.844. The 
inset shows the separation parameter cr as a function of e. The 
dashed line is the linear asymptote a = Ce The range where 
the asymptote deviates less than 10% from a is between the 
diamonds that are at ip = 0.676 and (p = 0.714. 



As mentioned in section II., the separation parameter 
cr(e) can be calculated from J-q[fk] and its derivatives at 
(fie- The result is given in the inset of Figure |5] The 



linear term in Eq. (17 1 describes a{e) for —0.030 < e < 
0.025 with an accuracy better than 10%. This range is 
similar to that for the corresponding result for d = 3 [25], 
and provides an estimate for the range of validity of the 
asymptotic expansions. The quality of the leading order 
result and the next-to-leading order contribution for 
(cf. Eq. (15l) is demonstrated in Figure 5 for the three 
q- values qi, q2 and ^3 (Figure 4a). It is interesting, that 
for 52 the leading asymptote describes fq{e) best and for 
an unexpectedly wide range. This arises because of a 
cancellation of higher order correction terms. For qi and 
93 the next-to-leading order has to be taken into account 
already for rather small e, as expected from the cr(£)- 
curve. The overall behavior also is quite similar to d = 3 
(cf. Figure 3 of Ref . ^22]). 



Now we turn to dynamical features. Figure 6 presents 
the normalized correlators (j)i{t) = 4>q. {t) for i = 1,2. The 
two-step relaxation process for e < becomes obvious for 
both correlators. Since /^^ < /^^ (cf. Figure 4a,b) the 
plateau heights for 02 (^) are below those for 4>i{t). Again, 
the t- and cj-dependence (the latter is not shown) is in 
qualitative agreement with the corresponding results in 
d = 3 (cf. Figures 4 and 6 of Ref. [l^)- In the following, 
we will apply the asymptotic expansions from Sect. II 
to the correlators in order to characterize the long-time 
dynamics in more detail. 
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FIG. 6: Coherent correlators 01 (t) (top) and 02 (t) (bottom) 
for different e = ilO""''^ Curves are labeled by n. 



Following Ref. [22] we have calculated typical time 
scales T^(e^O) and (e < 0) characterizing the first 
and second relaxation step. In the fluid, T^{q) marks 
the crossing of the plateau, (j)q{T^{q)) — and the a- 
relaxation time is defined by 4>q{T'^) — fql&- In the glass, 
T^{q) captures the approach to the long-time plateau, 
<l^q{r+) - fq = 1-001(/, - rq). Rcsults arc shown in Fig- 
ure [7| for q = Qi. The divergence of these relaxation times 
at e = is described by the asymptotic laws, Eqs. (20a) 
and ( |20b I. Since ip^^'^ has been determined, one can cal- 
culate the exponent parameter A. As a result we find 
j^d=2 0.7167 which implies a'^^'^ ^ 0.320, b'^=^ ^ 0.613 
and 7''=2 2.38. These values are close to A''^^ = 0.735, 
^d=3 ^ 0.312^^^=3 ^ 0.583 and -/'^=^ = 2.46 [22 . Us- 
ing in Eqs. (20a) and (20b I a'^^^ and 7''=^ leads to the 



asymptotes (solid lines) in Figure [t] 
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FIG. 7; Timescales for the f3 process and r' for the a 
process in the liquid for wave vector qi. The lines are the 
asymptotic power laws. A — 0.7167, a = 0.320, b — 0.613, 

= 0.009, + 2.38, 



S = j- = 1.56, c+ = 0.044, 



0.0173 (For the definitions of c± and c' see Ref. 



The microscopic time scale tg entering the critical 
power law, Eq. (18), can be deduced by plotting {(j)q{t) — 



fq)t°' versus log t for e > and e < (see Figure 8| 
The value at which a constant plateau is best reached d;' 

With hq^ ^ 0.337 and a'^=' 



both curves is 0.114= 
we get to = 0.034. 
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FIG. 8: Determination of to. {4>i{t) — ft)t°' is plotted over 
logt, for glass and liquid curves close to the transition. The 
position where the plateau is best reached by both curves is 
marked by arrows. The value found is 0.114. to is then given 



0.337 and a = 0.320. 



by (*^:i^)V- = 0.034, with hi 

hi 



The quality of the leading order result (Eq. ( [T9| ) and 
its next-to-leading order correction (second and third 
term in the curly bracket of Eq. ( 18 1) of the critical law 
is checked in Figure |9l The similar check for the von 
Schweidler law (Eq. 

cc 
J a 



\2l\ ) is done in Figure 10 Like for 
, ^ , the leading order has a large range of validity for 
q = 92- 
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FIG. 9: The critical laws. The leading asymptotes (dashed) 
describe the solution within 10 % up to the points marked by 
diamonds (41, 0.52) for {qi, 52)- The next to leading results 
(dotted) are within 10 % up to the circles (0.73, 0.17) for (qi, 
52)- Correction amplitudes are K(a) = —0.021, Ki — —1.007, 
K2 = -0.120. 



Note, that both time scales of the structural relaxation 
are given by the matching time to determined in Fig. |8j 
and by the separation parameter e. Thus, there remains 
no adjustable parameter in the test of the von Schweidler 
law in Fig. [T0| 




FIG. 10: The von Schweidler law. The leading asymptotes 
(dashed) describe the solution within 10 % up to the points 
marked by diamonds (0.29, 0.053) for (qi, q^). The next to 
leading results (dotted) are within 10 % up to the circles (0.37, 
0.47) for (gi, 52) {b = 0.613 and = 7.4510^). Correction 
amplitudes are K{b) = 0.496, Ki = -1.007, K2 = -0.120. 



describes the first scaling law regime. The t- and e- 
dependence of (f)q{t) = {(j}q{t) — fg)/{hqC) on the time 

scale tcr is given by g±{t) for |e| 0, independent on 
q For the two-dimensional system this property is 



demonstrated by Figure 11 
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FIG. 11: (Color onhne) Functions 4>q{t) = {(pqit) - 
fq)/{hqVC) for gi = 6.46, 92 = 10.06 and gs = 18.26 for 
£ = ±10""''^ for 3 values of n (solid lines). The scaling asymp- 
totes G{t)/\/CJ are shown as dashed lines. The region where 
the functions for different q collapse onto a master function 
(/3-region) increases with decreasing e. C — 2.08; the curves 
for £ > are shifted down by 1. 



The critical law and the von Schweidler law are the 
short and long time expansion of the so-called /3-master 
function g±{t = t/ta) for 0, respectively. g±{i) 



We clearly observe that the curves for qi, q2 and 
^3 collapse onto a master function with increasing n, 
i.e. for e ^ 0. Because of the connection between the 



q-dependences of the correction amplitudes in Eqs. (18 1 
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and (21 1, an ordering scheme exists for the functions 
(jyqit) in Fig. [ill Their vertical order before and after 
crossing the plateau needs to coincide; this is obeyed in 
Fig. [11] 

The second scaling-law regime (for e < 0) is defined 
by the rescaled time t = t/t'^, where t'^ often is called 
a-relaxation time and denoted by r. For i — 0(1), the 
t-and e-dependence of (j)q is given by the a-master func- 
tion (l)q{t) p" . The validity of this second scaling-law 



(Eq. (22l) is presented in Figure 12 Approaching e = 



from below a collapse onto a g-dcpcndcnt master function 
dn occurs. 




FIG. 12: (Color online) Correlators (t>q(t) for different 
£ = —10^"^^, n=4,6,8,10 as function of rescaled time 
t — t/t'^ (solid lines). The thick solid line is the a- 
master function </>q(t). The dotted lines are the short time 
parts of leading-plus-next-to-leading approximatio n fo r the 
Qf-correlators (f>q{i) + hqBiai^^ according to Eq. (221 with 
Si = 0.5/(r(l - b)r(l -h 6) - A) = 0.374. Light (magenta) 
curves give the Kohlrausch laws fitted to the 0q(t) in the range 



in Fig. 12 far exceeds the one of the /3-scaling law tested 
in Fig. 11. This originates from the dependence of the 
leading corrections on the separation parameter e. While 
the corrections to the /3-process are smaller by a factor 
only, the relative corrections to the a-superposition 
principle start out in order e. For example at (72, the a 
master function (j)q^ (i) describes 68% of the decay of the 
final relaxation better than on a 5% error level at the 
separation e — —0.01 (see the circle in Fig. 12 1, while 
in Fig. [11] for the test of the /3-scaling law smaller e are 
required. 

The shape of the a-relaxation process, viz. its master 
functions (pqit), often is described by a Kohlrausch law 

^,(i)«A, exp{^(^pj "} (31) 

where a possible dependence of the parameters on 
wavevector q is taken into account. The von Schweidler 



expansion of the a-process, Eq. (21 1, immediately shows 
that the <^g(t) exhibit stretching, viz. do not decay via 
simple exponential relaxation, but result from a broad 
distribution of relaxation times. Numerical solutions of 
the MCT equations in rf = 3 have shown that, e.g. for 
hard spheres, the Kohlrausch law provides a good overall 
fit to the master functions for all wavevectors albeit with 
noticeable deviations at short rescaled times, and with 
q-dependent parameters [33j . Interestingly least square 



fits of Eq. (31 1 to the (final part of the) a-master func- 
tions for t > lO"'^'^®, yield remarkably close agreement 
of the d = 2 curves with Kohlrausch laws. This is shown 
for two wavevectors in Fig. 12, and can be learned from 
the fit parameters presented in Fig. 13. 




0.4 < 



logjQ(t) = [—3.86, 2.14]; the parameters are given in Fig. 13 



For the test of the superposition principle, the a- 
relaxation time was computed using the power-law 
Eq. (20b I and to from Fig. 8. At qi, the a-relaxation time 



obviously deviates early from the asymptotic power law, 
since there are intersections of the rescaled correlators. 
Nevertheless, the range of validity of the a-scaling law 



FIG. 13 : Kohlrausch fit parameters of a least square fit of 



Eq. (31 1 to the MCT a-master functions 4>q(t) on logarithmic 
time axis in the interval logiQ(t) = [—3.86,2.14]. Aq (circles) 
is quite close to fq (thick dashed). The Kohlrausch exponent 
l3q (diamonds) converges to the von Schweidler exponent b = 
0.613 (thick line) for high q values. The times fq are shown 
by a dotted line. 
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The amplitude Ag of the Kohlrausch law closely fol- 
lows the critical NEP, which gives the (true) amplitude 
of the a-process, and the Kohlrausch stretching expo- 
nent f3q varies little with wavevectors and is quite close 
to its large-g limit given by the von Schweidler exponent 
b. Within MCT, the Kohlrausch law can be derived as 
limiting law when an increasing number of correlators, 
i. e. correlators with a large range of q- values, con- 
tributes to the memory kernels |34j . This arises for large 
wavevectors, then (3q = h holds, and the Kohlrausch- 
relaxation time depends on wavevector as fq ~ q~^/^ . In 
the opposite limit, when only one correlator contributes 
to the memory function, the schematic model is ob- 
tained, where the a-process is exponential [2]. The lat- 
ter description obviously best applies to the relaxation 
of the correlator (j)q-^ (t) at the position qi of the pri- 
mary peak in the structure factor. It corresponds to 
the motion of particles connected to their average sep- 
aration, which in MCT is predominantly coupled back to 
itself. Guided by these two limits, we speculate that the 
Kohlrausch law provides rather good fits to the a-process 
in d — 2 because the spread of Kohlrausch exponents be- 
tween 6 = 0. 61 < Pq < 1 is smaller in d = 2 than it is in 
d = 3. 




FIG. 14: (Color online) Scaled mean squared displacement 
Ar^(f) = 5r^{t)/{2d) in units of for disks in d = 2 (black 
solid) and spheres in d = 3 (green dashed) for different e = 
±10""/^ as labeled. The thick lines are the critical curves. 



Finally, we have calculated the mean squared displace- 
ment Ar^{t) = Sr'^{t)/{2d), weighted with l/(2d). Fig- 
ure 14 presents the results obtained from the solution of 
Eq. (|24| for d = 2 and d = 3 [23] . For the hquid phase the 
increase of if leads to the formation of a plateau which 
has its origin in the cage effect, independent on d = 2 or 
d = 3. The dynamical behavior for t small is governed 

(s) 

by the short time diffusion constant Dq and that for 



i ^ cx) by D{ip) ~ {t'^{ip))~^ ~ {ipc — , the long time 
diffusivity. At ipc, i.e. for e = a transition occurs where 
the cage has an infinite life time such that D{lp) — for 
Lp > ^Pc- Consequently the particle becomes localized 
with finite localization length (cf. Eq. ([27|). It is in- 
teresting that rf=2 0.077 D and rf=3 o.075 D are 
almost the same, as can be seen from Figure 14 (as well 
as from Fig. |3|. 



Although our results were obtained for monodisperse 
hard discs we have made a comparison of the mean 
squared displacement (MSD) 5r'^{t) with corresponding 
experimental results of the binary system studied in Ref. 
[19j . The control parameter varied in the experimental 
system is the interaction parameter, or inverse tempera- 
ture, r. Interaction parameters larger than Tm = 60 are 
used, the value at which the corresponding monodisperse 
experimental system forms a crystal. Because we aim 
for a qualitative test of our MCT results only, like the 
existence of finite localization length at the transition, 
we concentrate on the e-insensitve parts of the MSD's. 
Thus, we arbitrarily choose £ « for the experimental 
curve at the highest F available, F = 592. The result of 



this comparison is shown in Figure 15 



In the fit a global length- and timescale was obtained 
by fitting the F — 592 data. The length- and time scales 
found at this F are then used in fits to data for lower F. 
Quite good matching to the data is achieved as shown 
in Figure [15] for the example of F = 474. Hence the 
scaling factors for time- and lengthscales can be taken 
as constants and independent of parameter F. With the 
scales set, the only fitting parameter left is ^p, which is ad- 
justed by eye. We find the fits obtained remarkably good, 
considering the simplification to map the binary experi- 
mental system onto monodisperse hard discs. From the 
fitted length scale of the single particle MSD of the ma- 
jority particles (viz. the big ones) a reasonable value for 
the effective particle size follows. The position of the pri- 
mary peak in the hard-disc structure factor closely agrees 
with the position of the experimentally obtained (partial) 
structure factor of the big particles; see the upper inset 
of Fig. [15] We thus conclude that MCT correctly cap- 
tures the ratio of localization length to average particle 
distance. 



13 



10 




E 

zS. 






5- 


-(— • 
i_ 




to 

V 


/// 3- 


1 


- / // 2- 
/V/ 1- 



r=592 




0.3 



0.4 



r=474 



b+s^ 



-0.2 



200 



log(t) 



FIG. 15: Comparison of experimental data (dotted) from |35j 
with theoretical curves (solid). Since the experimental system 
is a binary mixture three MSD curves are measured. Curve b 
is measured considering only big particles, s only small parti- 
cles and b-fs takes all particles into account without discern- 
ing big or small. The top panel is at F = 592, e — —10"'', the 
bottom one at F = 474, e = -10"^. The data at F = 592 can 
in the given region also be fitted by a critical curve and could 
thus correspond to a state in the glass. At F = 592, the unit 
of time is fitted as — l/470s, and the unit of length is 
fitted as (420,700,500) (^m)^ for the big (curve b), the smaU 
(curve s), and all particles (curve b-l-s). This leads to a short 
time diffusion coefficient _Do ~ 0.9e~^''m'^/s and an effective 
diameter of -D^g^ ~ 20.5/im for the big particles. The lower 
inset shows the fitted separation parameter e versus F. The 
upper inset shows a comparison of the (partial) structure fac- 
tor of the big particles at F = 592 with the one of the hard 
disk system at £ = 0; fortuitously, deviations only appear at 
larger q. 



V. SUMMARY AND CONCLUSIONS 

In a first step we have determined the explicit depen- 
dence of the MCT functions T(^<^}^ and T^^\(^k, '^^k\ '^'^ 
the spatial dimensionality d. This has also been done 
in the hydrodynamic limit q ^ 0, which generalizes the 
well-known result for d = 3 [TJ [53] to arbitrary di- 
mensions. 



The major motivation of the present contribution has 
been the investigation of the existence of a dynamic glass 
transition in two-dimensional systems. As model sys- 
tem we have chosen monodisperse hard discs which might 
be a reasonable approximation for a system of polydis- 
perse hard discs, at least on a qualitative level. Taking 
the static input quantity Sq from a modified, hypernet- 
ted chain approximation we have found that an ergodic- 
nonergodic transition occurs at a critical packing fraction 
j^d=2 ^ 0.697. At this critical density, both, the fluc- 
tuations of the collective and self density simultaneously 
freeze into a glassy state, as for d = 3. That ^p^"^ is about 
35% above ^p'^'^ might be not surprising, since the pack- 
ing fraction (/JtrTang — tt/a/TZ ^ 0.9069 of the triangular 
lattice is also larger then (p^cif ' value for the hexag- 
onal closed packed lattice, by about 20%. The difference 
between ^f^r'^ and (/^^^'^ becomes even more obvious when 
scaling is done with the random close packing values 
(^,^=2 ^ 0.82 and ^^'^=^ ^ 0.64 {ipJip,,^Y=^ ^ 0.85 

and ((^c/'/'rcp)''"'^ — 0.81 deviate by not more than 5%. 
Consequently, ip'^^^ might be a reasonable scale for ip'j:. So 
far MCT applied to d = 2 provides an explanation for the 
dynamic glass transition observed in Refs. [TTJ [121 US] ■ 
Since MCT overestimates the glass transition it is not a 
surprise that « 0.80 is above (/^^^^ = 0.697, quite 
similar to d = 3 [3 . 

For all the investigated quantities and properties con- 
nected to the local 'cage motion' on an intermediate time 
window, like the critical nonergodicity parameter, critical 
amplitudes, the exponent parameter, the amplitudes Kq 
and Kq of the next-to-leading order corrections, etc. we 
have found a weak dependence on dimensionality only, 
comparing d = 2 and d — 3. The largest change be- 
tween d = 2 and d = 3 occurs for the coherent NEP at 
higher wavevectors. For the two-dimensional case, the 
q-dependence is more abrupt. Otherwise, the similar- 
ity also holds for the localization length which differ 
by less than 3%. Accordingly, the Lindemann criterium 
applied to the melting of the glass phase is almost d- 
independent, at least for d = 2 and d = 3. For longer 
times, we found the '(time-temperature-) superposition 
principle' of the a-process; the correlators collapse onto 
a non-exponential master function. Interestingly, we ob- 
served for hard discs that the functional form of the a- 
relaxation closely resembles a Kohlrausch law. This holds 
better in d = 2 than for hard spheres in d = 3. While the 
prediction of a superposition principle for the final decay 
is guaranteed by the general structure of the MCT equa- 
tions, the shape of the relaxation process provides infor- 
mation on the local particle rearrangements. Apparently, 
in the lower dimension the memory kernels arise from 
a large number of contributions so that the Kohlrausch 
law as a limiting law of large numbers provides a bet- 
ter approximation to the cooperative motion during the 
a-decay in d = 2 than in d = 3. There is another conclu- 
sion we can draw from our results. On a qualitative level, 
we have found consistency with the observations made in 
Refs. [Ul [TSJ [T7] . Particularly, the stretching found for 
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the monodisperse Lennard- Jones system [TB] and the bi- 
nary mixture of soft discs [17] as well as the two-step 
relaxation process |T7] can be described by MCT in two 
dimensions. 

Of course, more quantitative comparisons are neces- 
sary. MCT should also be worked out for binary hard 
discs without and with pair interactions. This will al- 
low to compare MCT results with the experimental ones 
[m in detail. A first attempt has been done concerning 
the mean squared displacement. A more or less satisfac- 
tory agreement has been found over about four decades 
in time (see Figure 15). Whether the systematic dis- 
crepancies between the theoretical and the experimental 
result can be attributed to the different two model sys- 



tems, i.e. on one side the monodisperse hard discs and 
on the other the binary hard discs with dipolar repulsion 
at rather low densities, is one of the open questions we 
intend to study in the future. 
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